Instability of QCD ghost dark energy model 
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We investigate the instability of the ghost dark energy model against perturba- 
tions in different cases. To this goal we use the squared sound speed v 2 s whose sign 
determines the stability of the model. When v 2 < the model is unstable against 
perturbation. At first we discuss the noninteracting ghost dark energy model in a 
flat FRW universe and find out that such a model is unstable due to the negativity 
of the v"l in all epoches. The interacting ghost dark energy model in both flat and 
non-flat universe are studied in the next parts and in both cases we find that the 
squared sound speed of ghost dark energy is always negative. This implies that the 
perfect fluid for ghost dark energy is classically unstable against perturbations. In 
both flat and non flat cases we find that the instability of the model increases with 
increasing the value of the interacting coupling parameter. 

I. INTRODUCTION 



X 



Nowadays, there are enough observational evidences which indicate that our universe is 
currently experiencing a phase of acceleration [l-5]. The unknown cause of this unexpected 
acceleration called "dark energy" (DE) whose nature is still a matter of much doubt. Trying 
to explain observed acceleration of the universe expansion two main approaches have been 
followed in the literatures. The first is based on the DE proposal which assumes that there 
exist an exotic form of energy with negative pressure which is responsible for the acceleration 
of the universe expansion. Famous alternatives in this approach including cosmological 
constant A, scalar field models of DE such as quintessence fj, |7|, K-essence |8l-llC|. phantom 



fields 



llj , tachyon 12[ , holographic DE [13| and agegraphic DE [14| and so on (for a recent 
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review on DE models see 



15| and references therein). The second approach for explanation 



of the acceleration expansion is based on the modification of the gravity theory. In this 
approach the DE problem is considered as a shortcoming of the Einstein's gravity and try 
to modify the standard model of cosmology in such a way that the phase of acceleration is 
reproduced without including any new kind of energy. Examples of these models are f(R) 
gravity [l6,[l3| and braneworld scenarios 18, 19]. 



Seeking a solution to U(l) problem , the so-called Veneziano ghost has been proposed in 
ihe low energy effective QCD where they are completely decoupled from the physical sector 
20 ^23]. The ghosts make no contribution in the flat Minkowski space, but make a small 



energy density contribution to the vacuum energy due to the off-set of the cancelation of 
their contribution in curved space or time-dependent background. This contribution to the 
vacuum energy density can be considered a possible candidate for the origin of the cosmo- 
logical constant 24j. In a dynamic background or a spacetime with non-trivial topology 
the ghost field contributes to the vacuum energy proportional to Aq CD H, where H is the 
Hubble parameter and Aq CD is QCD mass scale. With A QC d ~ lOOMeV and H ~ 10~ 33 eV 
Aq CD H gives the right order of magnitude ~ (3 x 10~ 3 eV) A for the observed DE density 



24j . This remarkable coincidence imp lies that the ghost dark energy (GDE) model is free 



from the fine tuning problem 24], |25| . It was shown that this vacuum energy density can 



27]- 



play the role of DE in the evolution of the universe |26l . 

Every new model of DE represents new features and consequences which should be ex- 
plored carefully. One way to test the viability of a new DE model is to explore its stability 
against perturbations. To investigate the stability, a key quantity is the squared speed of 
sound v 2 = dp/ dp [28|. The sign of v 2 plays a crucial role in determining the stability of the 
background evolution. If v 2 s < 0, it means the classical instability of a given perturbation. 
This issue has already been investigated for some DE models. It was shown that chaplygin 
gas and tachyon DE have positive squared speeds of sound with, v 2 = —w, and thus they 
are supposed to be stable against small perturbations [29|, |3Qj. However, the perfect fluid of 
holographic DE with future event horizon is classically unstable because its squared speed is 



always negative [31] . Also in 



321 ] . it is shown that the agegraphic model of DE have a neg- 



ative sound speed squared in flat, non-flat and also in the presence of interaction indicating 
the instability of this model against perturbations. 

DE and dark matter (DM) are usually considered as two distinct dark components of 
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the universe due to their unknown nature and very different origin of them. Note that DM 
has gravity effect and plays a crucial role in explanation of the galaxy rotation curve, while 
DE has anti-gravity feature which pushes the universe and makes it expansion accelerated. 
However, in recent years several signals have been detected, implying a small interaction 
between DE and DM is possible. As an instance, observational evidences provided by the 
galaxy cluster Abell A586 supports the interaction between DE and DM 33J. Beside that, 
lately there arose an interest in the non-flat version models of DE in literature. This new 
enthusiasm originate from some observations which challenge the idea of the flat universe. 
For example evidences from CMB and also supernova measurements of the cubic correction 



to the luminosity distance favor a positively curved universe 34 



3511 . In addition, some 



exact analysis of the WMAP data reveals the possibility of a closed universe 36 1. 

All above reasons motivate us to study the stability of interacting GDE model in a nonflat 
universe. In this paper, we would like to generalize the approach presented in [28|, |3l|, |32| to 



GDE model in a universe with spacial curvature in the presence of interaction between the 
dark matter and DE. Various aspects of GDE have recently investigated. A thermodynamical 
description of GDE is discussed in 37j. Tachyon and quintessence reconstruction of GDE 



model were studied in 1 381 and 

n 

Brans-Dicke theory [401 ] . 



respectively. The study has also been extended to 



This paper is organized as follows. In the next section, we review the GDE model in both 
flat and nonflat universe. In section IHIl we explore the stability of the GDE model in all 
discussed cases of section II. We summarize our results in section [TV] 

II. A BRIEF REVIEW ON GHOST DARK ENERGY 
A. Noninteracting ghost dark energy 



Let us at first review the noninteracting GDE in a flat FRW background filled with a 
matter component and GDE. We follow the method of 2^]. Dynamic of such a universe is 
determined by the Friedmann equation 

1 



H 1 



3M p 2 



(Pm + Pd) 



(1) 



where p m is the energy density of pressureless DM and po is the GDE density. According 
to standard cosmology one can define the fractional density of different energy components 



4 



of the universe as 

= — , = —, (2) 

Per Per 

where the critical energy density is p cr = 3H 2 M 2 . Thus, the Friedmann equation can be 
rewritten as 



n m + n D = 1. (3) 



Here, we consider the GDE which its energy density can be written as 24] 



p D = aH, (4) 

where a is a constant of order Aq CD and Aqcd is QCD mass scale. Based on the definition 
of fin and using (j4]), we get 

_ Pd _ a 

Energy conservation equations for different components read 

p m + 3Hp m = 0, (6) 
p D + 3Hp D (l + w D ) = 0. (7) 

Taking the time derivative of relation (j3J), and using the Friedmann equation (TjQ) as well as 
the continuity equation (J7J) we find 

^ = ^ = -^H[l + n D w D ]. (8) 
Pd ti I 

Substituting this relation in continuity equation ((7j) we obtain the equation of state (EoS) 
parameter of GDE, namely 

Evolution of wd versus Qd is shown in Fig HI One finds that at the early time where Qd <C 1 
we have wd = — 1/2, while at the late time where Qd — > 1 the GDE mimics a cosmological 
constant, namely wo = —1. Another interesting quantity to be calculated is the deceleration 
parameter 

i H 1 3 Q D 

q = - l -jp = 2-2j2^r (10) 

From above relation it is clear that at the late time where Qd 1, q = —1. At the early 
universe where — > we have q = 1/2, which corresponds to the matter dominated epoch. 
Evolution of q versus Qd is plotted in FigfJJ Two important points should be emphasized 
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FIG. 1: Evolution of wd versus £Id for noninteracting GDE. 
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FIG. 2: Evolution of deceleration parameter against tip for noninteracting GDE. 



for noninteracting GDE model. First, in this model there exist no free parameter. Second, 
wd cannot cross the phantom line. The evolution of GDE is governed by 



27] 



Q! D = 3Q 



D 



2 - tt D/ 

where the prime denotes the derivative with respect to x — In a. 



(11) 



B. Interacting ghost dark energy 

DM and DE are usually considered separately in the modern cosmology. However, there 
exist some observational evidences which support an interaction between them 4lJ-Thus, it 
seems meaningful to study the interacting version of GDE. To this goal we start with energy 
conservation equations 



Pm + 3Hp m = Q, 



(12) 



6 



p D + 3Hp D (l + w D ) 



-Q. 



(13) 



In these equations the form of Q has not been chosen exactly yet, however, we know that 
Q should be small and positive. A large and positive choice of Q does not lead to late 
time acceleration while a large and negative choice of Q will lead to an early domination of 
DE component which does not let the large structures to be formed in the universe. One 
should note that a choice of negative sign for Q implies decaying of DE to the cold DM. In 
this paper we choose the interaction term as Q = 3H 2 (po + pm) which can be written as 



Q = 3H 2 p D (l + ^[ 



In the case of interacting, the EoS parameter is obtained as 



23| 



1 



w D 



1 + 



2b 2 



2 — V ' 

Setting b = 0, W£> for the noninteracting case is retrieved. In the late time where — >■ 1, 
the EoS parameter of interacting GDE necessary crosses the phantom line, namely, wd = 
— (1 + 2b 2 ) < — 1 independent of the value of coupling constant b 2 . Also the deceleration 



parameter can be calculated as 27 ] 

_ i 3 n D 

q ~ 2 ~ 2 (2 - Q D ) 
The equation of motion of Qe> can be obtained as 



2b 2 

1 + n~ D 



D 



din a 



2 D 



n 



D 



(2 - Q D ) 



2tf_ 



(15) 



(16) 



C. Interacting ghost dark energy in non-flat universe 



In the past decade several observational evidences have been observed in contrast to the 
flatness assumption of the universe. In the context of inflation which assumes a large rate 
of expansion at the beginning instants of the evolution of the universe, it is argued that 
the flatness is not a necessary consequence of inflation if the number of e-folding is not 



very large 



421 ] . Besides, the parameter Qf. represents the contribution to the total energy 



density from the spatial curvature and it is constrained as —0.0175 < Qk < 0.0085 with 95% 



confidence level by current observations 



431 ] . These motivate us to study the GDE model in 



the presence of curvature. Taking the curvature into account, the first Friedmann equation 
is written as 



H 2 + - 2 

a 2 



3M 2 



(Pm + Pd) ■ 



(17) 
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This equation can be written as 



VL m + VLd — 1 + 



(18) 



where = -Jj^ . We can obtain different parameters in the presence of the curvature term. 
Taking the time derivative of the Friedmann equation ffTTj) and using ffTB"]) . we find 



H 3 

= Q k - - [1 + Q k + Q D w D ) 



(19) 



Combining the above equation with Eq. ( !T3|) . after using Q = 3H 2 po(l + W^), a s we h as 



-j| we obtain the EoS parameter as 

1 



n k 2b 2 



(20) 



The deceleration parameter can be calculated by substituting Eqs. ( 1T91 and (120]) . We find 



3fi 



Q 



D 



2 2(2 - n D 
Also the equation of motion of GDE can be obtained as 



tt k 2& 2 , 



(21) 



27] 



dlna ~2 D V 3 2-^ D 



(22) 



It is worth mentioning that all relations we obtained in this subsection restor their respective 
expressions in the previous subsection when we set Q k = 0. 

III. INSTABILITY OF THE GHOST DARK ENERGY 



The main idea for investigating the stability of GDE model comes from the perturbation 
theory. Assuming a small perturbation in the background energy density, we would like 
to see if the perturbation grows or will collapse. In the linear perturbation theory, the 
perturbed energy density of the background can be written as 



p(t,x) = p(t) + 5p(t,x) 



(23) 



where p(t) is unperturbed background energy density. The energy conservation equation 
(V M T^ = 0) yields [28] 

5p = v 2 s V 2 6p(t,x), (24) 
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where v 2 s = jf- is the squared of the sound speed. Solutions of equation tfMty include two 
cases of interest. First when v 2 is positive Eq. (j2l|) becomes an ordinary wave equation 
whose its solutions would be oscillatory waves of the form 5p = 8poe~ loJt+lk ' x which indicates 
a propagation mode for the density perturbations. The second is when v 2 is negative. In this 
case the frequency of the oscillations becomes pure imaginary and the density perturbations 
will grow with time as 5p = 5poe ult+lk ' x . Thus the growing perturbation with time indicates 
a possible emergency of instabilities in the background. 

In what follows, we are trying to see if the GDE leads to instabilities in the FRW back- 
ground. To this end we will calculate the sound speed for both interacting and noninteracting 
GDE model. 



A. Instability of noninteracting ghost dark energy 

In this section we would like to obtain the sound speed in a background filled with matter 
and GDE in the absence of interaction. The squared sound speed can be written as 

«- 2 = T P < 25 > 

In order to obtain the sound speed in a background filled with barotropic fluids we rewrite 



the definition as 26j] 



, dP P p 

v 2 = — = - = ^w + w, 26 
dp p p 

where in the last step we have used P = wp. Taking time derivative of Eq.(j4]) and using 
(J7J), we obtain 

* = = =? (27) 

H p 3H(l + n D w D ) V ; 

Inserting fl9]), into above relation we get 

p _ 2-VL D 
~p ~ ~3H(l-n D ) 

Taking the time derivative of (Q yields 



(28) 



*" = -<A- (29) 



Replacing Eqs. (l28l) and ( 1291) into ( 1261) and also using (TTTT) we get 
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FIG. 3: Evolution of squared sound speed v 2 against Qd for noninteracting GDE model. 



where we also used 4 = ff-rp— . This result is the same as one presented in 26|. Having 
the Vg at hand we are ready to discuss about the stability of perturbations. It is clear 
that in the evolution history of the universe Qd can achieve values in the range [0, 1] which 
Qd = indicates beginning stages of the evolution of the universe while Qo — 1 is an 
extrapolation of the future of the universe. One can easily see from (130]) that v 2 is always 
negative and varies between [— 1,0]. Hence, Eq. ([231 has solution of the second category 
(Sp = 5p e ujt+ ' lk ' x ). This result indicates that due to the negativity of the squared sound 
speed every small perturbation can grow with time which leads to an instability in the 
universe. Thus we cannot expect a noninteracting GDE dominated universe in the future 
as the fate of the universe. The evolution of the v 2 versus Qd is shown in Fig. [3j 



B. Instability of interacting ghost dark energy 



A same steps as the pervious section can be followed to obtain the squared sound speed 



v s for the interacting case. Taking time derivative of Eq. (TBI) we have 



n 



W D 



D 



(2-n D ) 2 



Ab 2 

i + S r<n D -i) 



D 



(31) 



Also from (El) one finds 



P 
P 



3H(l + Q D w D ) 



Taking into account relation = H-^- as well as ( 1T6|) yields 



dt d In a 

3 



(2 -n D )\ + n D 



(32) 



(33) 
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FIG. 4: This figure shows evolution of squared sound speed v 2 versus for interacting GDE 
model. The solid line corresponds to b = 0.05, dashed line to b = 0.2 and dashed-dot line to 
b = 0.35. 



Replacing these relations in f|26|) and after a little algebra one obtains 



-2- + 2b 2 



(2-n D ) 



(34) 



which is the squared sound speed for interacting GDE fluid. The evolution of v 2 against Q D 
is plotted in Fig. (|4j) for different values of the coupling parameter b. The figure reveals that 
v 2 is always negative and thus, as the previous case, a background filled with the interacting 
GDE seems to be unstable against the perturbation. This implies that we cannot obtain a 
stable GDE dominated universe. One important point is the sensitivity of the instability to 
the coupling parameter b. The larger b, leads to more instability against perturbations . 



C. Instability of interacting ghost dark energy in non-flat universe 



Finally, we study the instability of interacting GDE model in a universe with spacial 
curvature. From Eq. ffl9]) we can obtain 

-2 



P 



P 3H 



1 



Taking the time derivative of Eq. ( 120]) . yields 



n 



w D 



D 



2-tt D 



1 



2-o, ^-T + £< 1 + *>) + X< 1 + *> 



(35) 



(36) 
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FIG. 5: This figure shows evolution of squared sound speed v 2 versus for interacting GDE 
model in a nonflat background. The solid line corresponds to b = 0.05, dashed to b = 0.2 and 
dashed-dot line to b = 0.35. 

Having the above relations at hand we are in a position to obtain the squared sound speed 
v 2 . Replacing Eqs. (135]) . f )36|) in fl26|) . after using f j22|) . one gets 

f>2 



2(1 - n D ) 2 i - n D 2b 2 3n D -A 



(37) 



(2-fi D ) 2 3(2-fi D ) 2 fi D (2-fi D ) 2 
Setting Qk = = b the above relation reduces to the flat noninteracting respective relation. 
Also the squared sound speed of the flat interacting case can be retrieved when Qk = 0. 

In order to obtain an insight on the stability issue of the interacting GDE in a nonflat 
FRW universe, we have to discuss on the sign of v 2 during the evolution of the universe. To 
this end we plot v 2 s versus Qp, where the value of Qd indicates different epoches of evolution. 
The result can be seen in Fig. |5] which clearly indicates an almost same behavior as the flat 
interacting case. For all epoches the squared sound speed, (v 2 ), is negative indicating the 
instability of the universe against perturbations in the GDE background. Once again we see 
the crucial dependence of the stability to the coupling parameter b. Increasing b will result 
more instability in the universe which is clearly seen in Fig. [5j As a result the universe 
filled with DM and GDE even in the presence of the curvature cannot lead to a stable GDE 
dominated universe. 



IV. SUMMARY AND DISCUSSION 

Every new DE model should at first explain the acceleration of the universe expansion 
whose EoS parameter satisfies wd < —1/3. If the model passed this test, then it should 
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be investigated for further features and consequences. Since we know from observation that 
our stable universe is experiencing a phase of acceleration due to the domination of DE 
(assuming the DE approach for explanation of the acceleration of the universe is correct), 
thus every new DE model should be capable to lead a stable DE dominated universe. Thus 
the investigation on the stability of a proposed DE model is well motivated. Among various 
DE model the so called QCD ghost DE model was recently proposed to explain the DE 
dominated universe within the framework of standard model of particle physics and general 
relativity. It was argued that the vacuum energy of the Veneziano QCD ghost field in a 
time-dependent background can play the role of DE with an energy density proportional 



to the Hubble parameter H 24j. The advantages of this new proposal compared to the 



previous DE models is that it is totally embedded in standard model so that one needs not 



to introduce any new parameter, new degree of freedom or to modify general relativity (26]. 

In this paper we have explored the stability of the GDE model against perturbations 
in the flat/nonflat background and with/without interaction. We used the squared sound 
speed (Vg = ^) as the main factor for studying the stability. If v"l is positive the GDE 
would be stable against perturbations. When v\ is negative we encounter the instability in 
the background spacetime. We have discussed several cases including whether there is or not 
an interaction between DM and GDE and whether there is or not a curvature term in the 
background metric. Interestingly enough, we found that the GDE model is always unstable 
against perturbations. As a result, the universe filled with DM and GDE component cannot 
lead to a stable GDE dominated universe. We also observed that the instability of the 
interacting GDE increases with increasing the interacting coupling parameter b. 
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